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ABSTRACT 

We present J, H and K$ photometry of the Orion Nebula Cluster obtained at 
the CTIO/Blanco 4 m telescope in Cerro Tololo with the ISPI imager. From 
the observations we have assembled a catalog of about ~ 7800 sources dis- 
tributed over an area of approximately 30' x 40', the largest of any survey deeper 
than 2MASS in this region. The catalog provides absolute coordinates accu- 
rate to about 0.15 arcseconds and 3a photometry in the 2MASS system down 
to J ~ 19.5 mag, H ~ 18.0 mag, K s ~ 18.5 mag, enough to detect planetary 
size objects 1 Myr old under Ay ~ 10 mag of extinction at the distance of the 
Orion Nebula. We present a preliminary analysis of the catalog, done comparing 
the (J-H,H-K S ) color-color diagram, the (H, J-H) and (K S ,H-K S ) color-magnitude 
diagrams and the JHKs luminosity functions of three regions at increasing pro- 
jected distance from the Trapezium. Sources in the inner region typically show 
IR colors compatible with reddened T Tauri stars, whereas the outer fields are 
dominated by field stars seen through an amount of extinction which decreases 
with the distance from the center. The color-magnitude diagrams make it possi- 
ble to clearly distinguish between the main ONC population, spread across the 
full field, and background sources. The luminosity functions of the inner re- 
gion, corrected for completeness, remain relatively flat in the sub-stellar regime 
regardless of the strategy adopted to remove background contamination. 

Subject headings: open clusters and associations: individual (Orion Nebula Cluster) 
- stars: pre-main-sequence — stars: luminosity function, mass function 
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Introduction 



The Orion Nebula hosts the richest cluster of young (r ~ 1 Myr) Pre-Main-Sequence 
(PMS) stars within 1 kpc of the Sun and th erefore represents an ideal labora tory to 



understand the process of star formation (see 



Muench et al. 



2008 



O'Dell 



2008 



, for recent 



reviews). The cluster, both rich (n 2000 members) and dense (about 2 x 10 4 sources per 
cubic parsec at its center), is dominated by a small number of massive OB stars mostly 
clustered in the Trapezium multiplet (9 l Ori). Their UV emission has created a blister 
HII region (Orion Nebula, M 42, NGC 197 6) whose ioni z ation front is still carving the 



underlying Orion Molecular Cloud (OMC-1, 



O'Dell 



2008 



O'Dell et al. 



2009|). About half 



of the young cluster members, surrounded by their original circumstellar disks, have been 
already exposed to the hard-UV radiation of the Trapezium stars, whereas the other half 
remain enshrouded within the OMC-1, together with newer active sites of star formation. 

Visible and near-IR data, both in imaging and spectr oscopy, are needed to characterize 



1993). These 



the main physical parameters of the cluster population (Hi llenbrand! 
observations, however, are hampered by the brightness and non-uniformity of the nebular 
background. To overcome these difficulties, we have exploited the unique combination 
of sensitivity and spatial resolution offered by the Hubble Space Telescope (HST) to 
obtain accurate photometry of the cluster, especially at substellar masses (HST Treasury 
program GO-10246). Our extensive HST survey has been complemented by ground-based 
observations, imaging the Orion Nebula from the (7-band to the Kg-band at La Silla and 
Cerro Tololo (CTIO) observatories. The ground-based observations, carried out in parallel 
on the same nights (but at a different epoch than the HST observations), complement 
the deep HST data which saturate at relatively low brightness levels. Their simultaneity 
makes the derived stellar colors largely immune to the uncertainties associated with source 
variability. 
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The ground-based survey at visible wavelengths, performed with the Wide Field 
C amera (WFI) a t the ESO/MPG 2.2 m telescope at La Silla, has been recently presented 



m 



Da Rio et al 



(120091 ) . In this paper we present the ground-based IR survey, performed 
at the CTIO/Blanco telescope with the ISPI imager. Its combination of sensitivity and 
field coverage (about 30' x 40') ideally complements the previous surveys of this region. In 
particular, it reaches fainter magnitudes than the wide field survey (30' x 45') by Hillenbrand 
et al. (1998), while it covers a much larger area than the deep surveys done with Keck 
(HX s -bands, 51 x 5(1, Hillenbrand & Carpenter, 2000), HST-NICMOS (Jff-bands, 2'.3 x 2% 
Luhman et al. 2000), UKIRT (JJJf-bands, 36 arcmin 2 , Lucas & Roche 2000), NTT (JHK S , 
5' x 5', Muench et al. 2002) and Gemini (JHKs-bands, 26 arcmin 2 , Lucas et al. 2005). All 



these deep surveys concentrate on a fi eld centered around t 



immediate surroundings in the case of 



Lucas k Roche! (|2000j). 



re Trapezium stars, or in its 



In Section [2] we discuss the observing strategy, while in Section [3] we describe the data 
reduction and photometric calibration. In Section 4 we present our completeness analysis. 
The resulting source catalog is presented in Section 5, whereas in Section 6 we discuss the 
color-color and color-magnitude diagrams and the luminosity function derived from our J, 
H, and K^-band photometry. 



2. ISPI observations 

The Infrared Side Port Imager (ISPI) is the facility infrared camera at the CTIO 
Blanco 4 m telescope. ISPI uses a 2k x 2k HgCdTe HAWAII-2 array, with reimaging optics 
providing a scale of 0.3 arcsec/pixel corresponding to a 10.25 x 10.25 arcmin field of view. 
Our target area, about 30' x 40', covers the field imaged with the HST. 

The observations were performed on the nights of 1 and 2 January 2005 (indicated 
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hereafter as Night A and Night B respectively) in the J, H and K$ filters using the Double 
Correlated Sampling readout mode. The main filter parameters are listed in Table [TJ The 
seeing was ~ Of! 7 in most of the K^-band images and occasionally worse on night B. 

The survey area was divided into eleven fields (Figured]). The central position of each 
field, together with the airmass of the observations, is listed in Table [21 Each field was 
observed with an ABBA pattern. After pointing the telescope to the center of a field (A 
position) a 5 point dithering pattern was executed offsetting the telescope approximately by 
±30" along the field diagonals. This first group of 5 dithered exposures was followed by a 
second group of sky frames (B position) taken on a nearby field free, or nearly free, of diffuse 
nebular emission to extend the field coverage and monitor the flat-field response close in 
time to our observations. The sky-source sequence was then repeated back, completing the 
ABBA cycle. This provides a total of 10 frames of 30 s exposure for the common part of 
each field, totalling 300 s integration. Due to the high background, we split the Kg-band 
30 s exposures in two consecutive 15 s exposures, maintaining the background well within 
the linear regime. 

Each dithering sequence was repeated with a 3 s integration time. These short frames 
were coadded into another image of 30 s total integration time, which we used to extract the 
photometry of sources saturated in the longer exposures. Due to an error in our observing 
procedure, the J-band observations of field 3 were obtained only with 3 s exposures. 

For absolute calibration we observed at various airmass the faint IR standard stars 
9108, 9118 and 9133 of 



Persson et al. 



(jl998l ). Unfortunately, on both nights the atmospheric 
transmission at IR wavelengths turned out to be unstable due to variations in the water 
vapor opacity, and we eventually derived the zero point calibration of our images through 
comparison with the extensive set of 2MASS data across our wide field described 
in Section [3721 
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of the detector, w. 
N. Van der BliekE 



3. Analysis 

The first step in the data reduction was the correction for the intrinsic non-linearity 
diich we performed by applying the correction curve kindly provided by 
We flat-fielded each image using dome flats and correct for residual 
color terms at low-spatial frequency using delta-flats derived from median averaged sky 
images. The dome flats were also used to flag bad pixels. The sky frames, properly filtered 
to remove spurious sources arising from the latent images of bright stars on the infrared 
detectors, were combined and subtracted from the science images. 

Since ISPI exhibits significant field distortion, each science image had to be 
geometrically corrected before being combined in a dithering group. This required a first 
pass of aperture photometry with DAOPHOT, in order to extract a source catalog with 
well measured centroids and photometry (to disentangle the case of multiple candidates 
in the search radius). These sources were then cross-identified with the 2MASS catalog 
for astrometric reference, building an image distortion map relating their position on 
the original images to the actual position provided by 2MASS. Due to the non-uniform 
distribution of sources, clustered at the center of the nebula, the accuracy of the distortion 
correction map varies across each field. Therefore, we combined the distortion maps 
relative to all images taken on the same night to derive, for each filter, a master distortion 
map. Each individual image was then warped using a 4th-order polynomial fit with 
coefficients derived from the corresponding distortion map. Images belonging to the same 
dithered groups were then combined per integer pixels into single 2500 x 2500 images. The 
DAOPHOT source extraction was then repeated on these geometrically corrected images to 
produce a new source catalog from which we derived the final astrometric solution of each 



1 now also available available on the ISPI web page: 



http:/ /www. ctio.noao.edu/instruments/ir_instruments/ispi/ 
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image by minimizing the residual shift and rotation with respect to the 2MASS positions. 
The average scatter between our coordinates and the corresponding 2MASS coordinates 
turns out to be about 0.15 arcsec, i.e., half of a ISPI pixel. The same procedure has been 
adopted on the 3 s images. 

In conclusion, for each field the final data consist of one combined image in each of the 
3 filters (J, H, K s ) and exposure times (300 s and 30 s), with the exception of fields 4, 5, 
and 6, which have double this number of images, and field 3, which was not observed with 
deep J images. 

3.1. Instrumental magnitudes 

The extraction of source photometry was performed on the final combined images, for 
both short and long exposure times. The Daophot FIND IDL procedure was used to extract 
an initial list of candidate sources, with rather loose detection thresholds and criteria to 
include the largest possible number of candidates. This initial list, counting more than 
60,000 candidates, has been visually inspected on the original images to clean up artifacts 
and to preliminarily classify real sources as either point-like or extended. The list was then 
reduced to 7563 sources, divided between 6305 point-like and 933 diffuse. 

For all sources we first derived an aperture magnitude by integrating the flux over a 5 
pixel (1.5 arcsec) radius circular aperture, taking sky annuli typically between 10 and 15 
pixels (3.0-4.5 arcsec); for 1416 bright sources, saturated in the 30 s images, we used the 
photometry extracted from the 3 s images. 

We then performed PSF photometry on all point-like sources. Due to the spatial 
variability of the PSF across the field of view of ISPI, present even after geometric 
correction, we divided each image in 9 parts (3x3 squares) allowing for some overlap 
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between adjacent sub-images. The samples of stars in each sub-image, still relatively 
rich, had a more homogeneous PSF allowing us to reliably derive PSF magnitudes. We 
subtracted the PSFs and looked at the residual images searching for faint companion 
stars previously undetected. We found in this way 325 "hidden companions" . Finally, we 
performed a second pass of PSF photometry, first by removing from the list of stars used to 
derive the PSF those having a faint companion, and then deriving also the photometry of 
the newly detected faint companions. 



3.2. Absolute calibration 

As mentioned in Section [2} the IR photometric quality of our nights turned out to be 
unsatisfactory due to hygrometric variability. We therefore calibrated our ISPI instrumental 
magnitudes directly to the 2MASS photometric system. The 2MASS catalog provides 
absolute photometry in the J (1.25/im), H (1.65/im), and K$ (2.17/im) bands to a limiting 
magnitude of 15.8, 15.1, and 14.3, respectively, with signal-to-noise ratio greater than 10. 

To estimate the color terms and zero points of each image, we first scaled the 
instrumental PSF magnitudes of point sources measured in the short 30 s images to the long 
300 s images, chosen as the reference time because of their source richness. To minimize 
systematic effects due do atmospheric variations, we did not simply add the nominal 
Amag = 2.5 factor but compared the instrumental magnitudes of well measured stars to 
derive for each field and filter a mean magnitude offset Amag between the short and long 
images. By homogenizing in this way the zero points of the short and long images, we add 
a small uncertainty to the photometric errors of the bright stars measured in the short 
exposures. 
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For each bandpass A we then evaluated the relations 

ma#2MAss(A) - ma#ispi(A) = ZP X + e Al ,A 2 • Prspi(Ai, A 2 ), 

where mag 2 MAss is the magnitudes in the 2MASS photometric system, magispi is the ISPI 
instrumental magnitude, the intercept ZP\ is the zero point in a given bandpass, the slope 
e\ 1} \ 2 is a color coefficient between the wavelengths Ai and A2 and cigpi is the corresponding 
observed ISPI color. We found the strongest statistical correlation in the combination: 

-^2MASS — -f^ISPI = ZPh + €jk ■ ( J — i^)ispi. 

For each field, we derived the linear fit parameters ZPjj and €jk using an iterative 
procedure to rejects spurious outliers and all sources with errors greater than 0.1 m in the 
2MASS catalog and 0.05 m in our input catalog. As shown in Tabled the color coefficient 
€jk varies slightly from field to field during each night, whereas the zeropoint ZP\ shows 
an increase towards the middle of the night (observation were scheduled with fields 4 and 
10 transiting at the meridian, in the two nights). 

For the J and K$ bands, we derived the zero-points and the color terms from a linear 
fit to the color relations 

{J — -^)2MASS = ZP JH + €jh ■ (J ~ -H") ISPI 5 

{H — K s ) 2 mass = ZP H k + £hk ■ {H — Kg)ispi- 

This allowed us to derive magnitudes calibrated in the 2MASS system for all sources with 
ISPI JHK instrumental photometry. In Figure [2] we show a comparison between the 
magnitudes reported in the 2MASS catalog and our calibrated magnitudes. The plots 
show the lack of systematic differences, besides the random errors that we attribute to 
photometric errors or stellar variability. 

For sources lacking one or two magnitudes (typically either J or both J and H) , we 
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assumed a linear relation of this type 

mag 2 MAssW = a x rnagispi(X) + (3 

deriving a and /3, for each bandpass and field, from the sample of calibrated sources. The 
parameter a provides a magnitude-dependent correction to the zero point j3, which is 
appropriate since fainter sources have, on average, redder colors. Using these coefficients we 
calibrated the magnitudes, with the relative uncertainties, of all remaining sources. 

Figure [3] shows the photometric errors plotted as a function of magnitude for all 
point-like sources. The multiple threads in the distribution of the photometric errors are 
due either to the superposition of deep and short exposures (bright stars, measured in short 
3 s exposures, turn out to have relatively larger errors) or to the sources measured, with 
higher photometric errors, in the crowded central part of the cluster. For comparison we 
also plot the average mag vs. Delta mag relation for the 2MASS sources falling in our field. 

4. Completeness 

To determine the completeness of our photometric catalog, we used an artificial star 
experiment. For each field and filter, we averaged the 9 PSFs created for photometric 
extraction (Section 13.11) into a common reference PSF valid for the full image. We scaled 
the reference PSF in steps of 0.1 mag covering the entire magnitude range measured across 
each field and injected each scaled PSF at a random position in the image, with the only 
caveat of leaving enough distance from the border to allow for a meaningful measure. 
We assumed as a detection criterion a photometric error less than 0.3 mag and less than 
0.5 mag difference between the injected and recovered magnitude, making also sure that 
the stars recovered by the code were the artificial ones and not previously known real stars. 
The completeness was calculated as the fraction of successful detection after 1000 iterations 
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of the process, for each magnitude bin. 



It turns out that the completeness (and the sensitivity) of our survey is significantly 
affected by the nebular brightness. We have therefore defined three concentric regions 
(Figure H]) at increasing distance from ^Ori-C. The limiting radii, respectively 6.7', 14.3' 
and 27. 2' (correspond i ng to projected distances of 0.81, 1.72 and 3.27 pc at a distance of 



414 pc, 



Menten et al. 



(120071 )) have been set in such a way that the three regions contain 
the same number of point-like sources. In order to avoid significant contamination from 
the M43 cluster, we neglect in our analysis a circle with a 5' radius centered on NU Ori 
(RA(2000.0)=05 h 35 m 31.37 s , DEC(2000.0)=-05°16'02"6, see Figure SJ). 

The results of our simulations are shown in Figure [5j For each filter, the completeness 
estimated in the outer region is slightly, but systematically, larger than that in the 
intermediate region, which in turn is much larger than in the bright inner region, the most 
heavily affected by crowding and diff use emission, where our s e nsitiv ity b egins to drop a t 



HKs ~ 14 mag. A comparison with 



Hillenbrand & Carpenter! (120001 ) and 



Muenchetal 



(120021 ). who surveyed the inner part of the cluster (5' x 5' and 6' x 6' fields, respectively) 
deriving HK photometry down to our photometric limit (Kg ~ 18 mag), shows that we 
detect as point sour ces more than 90% of thei r sourc es (648/706 in the if-band and 657/698 



in the K^-band fo r 



for 



Muench et al 



obtained by 



Hillenbrand k Carpenter! (12000h and 601/662(H) and 665/714(K s 



2002')'). Figure El where w e ov erplot the Luminosity Functions (LFs 



Hillenbrand &: Carpenter! (120001 ) and 



Muench et al. 



(120021 ) (dashed lines) to 



our observed LFs (gray area) constructed for their same regions, confirms that the missing 
sources are generally fainter than HKs — 14 mag (the missing bright sources, detected 
but not measured by ISPI due to saturation, are not relevant here). If we apply to our H 
and Kg LFs the appropriate complete ness corrections, estimate d by re peati ng; the artificia l 



star experiment on the same fields of 



Hillenbrand fc Carpenter! (120001 ) and 



Muench et al 
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( 120021 ). our LFs become consiste nt with theirs, with the exception of the secondary peak 



at K s — 15.5 mag discussed by 



Muench et al. 



(120021 ) . In Section 16.31 we will apply the 



completeness corrections to the measured counts, when we will compare the LFs at various 
distances from the cluster center. 



5. Results 

Our final photometric catalog contains 7759 sources. Having visually inspected all 
sources, we find that 6630 can be classified as point-like sources and 933 as diffuse objects. 
The other 174 sources, typically brighter than K s — 10, turn out to be saturated even 
in our short 3 s images. We included them in our catalog, adopting the photometry and 
coordinates reported by 2MASS. We also found that the 2MASS photometric catalog misses 
22 sources in the Trapezium region clearly visible in the 2MASS images an d saturated eve n 
i n our short exposures. We have added their JHK$ photometry using the 



Muench et al. 



(120021 ) photometry. 



In Table H] we summarize the number of sources measured in different combination of 
filters, exposure times, and morphological types. The source catalog, presented in Table EJ 
is available in the electronic edition of this paper. In Appendix |A] we illustrate the content 
of each field of the source catalog. 

Figure [7] shows a JHK$ color composite of our imaged field. The image, produced by 
the graphic staff of the Space Telescope Science Institute, has been artificially enhanced 
to reduce the original dynamic range and more clearly reveal the inner part of the cluster. 
Figures [8]-19 provide 12 finding charts, starting from the Trapezium region (Figure E|) and 
continuing with the 11 fields (Figure 9 to 19). The sources are coded with circles having 
sizes proportional to their brightness and colors scaled from blue to red according to their 



observed J-K$ colors. 



6. Discussion 

In this paper we present a preliminary analysis of the data, concentrating on the 
point-like sources of the Orion Nebula Cluster having the full set of JHK$ photometry with 
magnitude errors Am < 0.3, with the exception of sources seen only in Field 3 (thus lacking 
the deep J exposure) and within 5' from NU Ori, i.e., lying in the field of the satellite 
nebula M43, about 15' northeast of the Trapezium. Our sub-sample is therefore composed 
by 4103 sources. 



6.1. Color-color diagrams 



In Figure 120a we show the near-IR J-H versus H-K$ color-color diagram. Following 



LeeetaL 



(120051 ). we overplot the intrinsic colors of the main-sequence stars and giant branch, 



togeth er with the 



disks (IMeyer et al. 



ocus of classical T Tauri stars (CTTSs) surrounded by circumstellar 
1997( 1. given by the relation (J-H) - 0.630(ff- Kg) - 0.497 = 



in the 2MASS photo metric system. We also plot the interstellar reddening vector of 



Cardelli et al. 



(119891 ). given by the relation E(J - H)/E(H - Kg) = 1.83 in the 2MASS 
photometric system. Accounting for extinction, the locus of the CTTSs lies between the 
two parallel dotted lines defined by the equations (J — H) — 1.83(ff — K s ) + 0.098 = and 
(J-H)- 1.83(H- Kg) + 0.50 = 0. 

For 45.2% (1853/4103) of the sources, the near-IR colors are consistent with those 
of stellar dwarfs (luminosity class V), giants (luminosity class III) or CTTSs reddened by 
various amounts of foreground extinction. For 11.4% (466/4103) of the sources the near-IR 
excess is compatible only with the reddened CTTSs, whereas only 0.8% (35/4103) of the 



-14- 



sources show near-IR excess significantly greater than those of CTTSs, indicative of strong 
K^-band excess. 

To assess how the observed reddening and color excess vary with the distance from the 
center of the cluster, we consider the Hess diagrams for the three radial regions. The inner 
field (Figure I2"0"b) shows a peak well detached from the main sequence, with iso-density 
contours elongated in a direction parallel to the CTTSs locus, together with a long tail of 
stars with high reddening. The Hess diagram for the intermediate region (Figure [20b) shows 
a main peak much closer to the locus of dwarf stars, together with a secondary maximum at 
higher extinction. The intermediate contour levels (e.g. the green area) are now elongated 
along the range of reddened main sequence stars. Finally, the Hess diagram for the outer 
region (Figure l2~0~bl) shows again a strong peak compatible with dwarf stars affected by a 
small amount of extinction, whereas the low contour levels appear less pronounced than in 
the two other regions and well within the limits of reddened main sequence or giant stars. 

Overall, the three color-color diagrams can be interpreted as follows: 1) the inner 
Trapezium region is mostly populated by young stars with near-infrared excesses typical of 
CTTSs. They have generally relatively modest amounts of extinction, which is compatible 
with the fact that the ionizing radiation from the Trapezium stars has cleared the molecular 
cloud and exposed them to our view; 2) the intermediate region contains a relatively larger 
number of heavily reddened sources. They are seen through denser parts of the OMC, not 
yet reached by the expansion of the ionized cavity; 3) the outer regions mostly contain 
background objects seen trough the low opacity layers at the ed ges of the main OMC r idge, 



which crosses the ONC approximately in north-south direction ([Goldsmith et al. 



1993). 
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6.2. Color-magnitude diagrams 
In Figure [21] we present the (if, J-H) and (Kg, H-Ks) color-magnitude diagrams still 



for th e same sample. Together with our data, we plot the 1 Myr isochrone of 



Chabrier et al. 



(12000 



CBAOO), which covers the range 0.001 M < M < 0.1 M & and is tailored to model 



substellar and planetar y mass obje c ts in the limit of high opacity photospheres. We also plot 



Siess et al. 



(12000 



the 1 Myr isochrone of 
for stellar sources. Both isochrones have 
using the transformations of 



sieOO) ranging from 0.1 M to 7 M , well-suited 
jeen converted to the 2MASS photometric system 



Carpenter! (120011 ). The offset between the two isochroness 



can be regarded as a graphic representation of the uncertainties in theoretical models. An 
Ay = 10 reddening vector is also indicated, as well as the reddening lines starting from 
the isochrones in correspondence of a 1 M® object (J=10.50, H=9.88 and X5=9.76 at 
414 pc in the lMyr SieOO isochrone), a M=0.075 M Q object at the hydrogen-burning limit 
mass (J=14.29, H=13.86 and Kg=13.56 at 414 pc in the lMyr CBAH00 isochrone) and a 
M = 0.012 M & substellar object at the deuterium-burning limit mass (J=17.33, ff=16.90 
and Ks=16.48 at 414 pc in the 1 Myr CBAHOO isochrone). For this distance and isochrone, 
we are therefore sensitive to stellar sources down to Ay ~ 60, brown dwarfs down to 
Ay ~ 40, and planetary mass objects down to Ay ~ 15. 

Almost all data points lie to the right of the isochrones and are therefore compatible 
with reddened 1 Myr old sources (as well as heavily reddened field objects). If we discard 
the blue sources with J — H < 0.3, typically affected by large measurement errors, we find 
that the reddening lines identify 2246 sources in the region of reddened stellar photospheres, 
1298 sources in the region of reddened brown dwarfs and 142 sources in the region of 
reddened planetary mass objects. In the (Kg, H-Kg) color-magnitude diagrams these values 
are 2134, 1145, and 421, respectively. 

Splitting the cluster in the same three regions of Figure HI we obtain the color- 
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magnitude diagrams presented in Figure [22j From top to bottom, Figure [22] shows the 
(H, J-H) (left column) and (K$, H-K$) color-magnitude diagrams for the inner, median 
and outer regions, respectively. Consistent with what was found analyzing the color-color 
diagrams, the inner and outer regions appear dominated by two different populations. 
The top diagrams clearly reveal the location of the ONC: an ensemble of stellar sources 
with modest amounts of foreground extinction. The bottom diagrams are dominated by 
fainter sources (of course, for background stars the 1 Myr isochrone at 414 pc does not 
apply), whereas the intermediate region contains a mix of both populations. A careful 
look at the diagrams show that the peak corresponding to the ONC drifts down (i.e., to 
lower masses) and left (to lower extinction values) moving from the inner to the outer 
regions. The shift to lower masses represents an overabundance of higher mass sources at 
the center of the cluster, i n agr eement with the evidence for mass segregation found by 



Hillenbrand fc Hartmannl (119981 ). The redder color of the cluster members in the central 



field may be due to the fact that they are, on average, either more embedded in the OMC or 
subjec t to higher circumstellar extinction fr om their heavily warped or flared circumstellar 



disks ( IRobberto. Beckwith fc Panagia 



2002j), or a combination of the two. 



6.3. Luminosity Functions 

We conclude our analysis showing the J, H and K^-band luminosity functions (LFs) 
for all point-like sources. Once again, we remove stars less than 5' projected distance 
from NU Ori and divide the remaining sample in three groups according to the distance 
criteria adopted in the previous sections. Rather than assuming a certain completeness 
limit, we fully apply our completeness corrections, derived in Section 0], to the source counts 
measured in the various regions. 

The histograms for the inner region (Figure [231 thick solid line) show a broad peak at 
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J ~ 13 , H ~ 12 a n d K,^ ~ 12, corresponding to an M6 star (0.175 M Q ) at 1 Myr, according 



to the 



Siessetal 



(120001 ) model with zero reddening. The peak shows an abrupt drop at 



about JH ~ 14 and K$ — 13, approximately by a factor ~ 2, w hich is fully consistent wit h 



the pr evious near-infrared studies of the ONC, e.g. Figure 5 of 



Hillenbrand fc Carpenter 



((20001). On the other hand, the luminosity functions corrected for completeness remains 
remarkably flat across the substellar region down to our sensitivity limit. This is not what 
has been found by other authors who concentrated on the innermost area. We remark that 
our completeness correction has the largest effect in the inner region, the one more affected 
by crowding and by the brightness of the nebular background. On the other hand, we have 



seen that our comp 



Muench et al. 



eteness correction fails to reproduce the K$ — 14 secondary peak of 



(120021 ). and for this reason can hardly be suspected of providing an artificial 



overcorrection in the central area. 

Moving from the inner to the outer parts of our survey area, the main peak in the LFs 
becomes less prominent in the median region and eventually disappears in the outer regions 
giving to the LFs an increasing monotonic slope. If we compare our Ks'-ban d LF for the 



outer cluster with the .fT-band LF of the off-cluster control field presented by 



Muench et al. 



(120021 ) (Figure 3a of their paper), which shows no stars brighter than K ~ 12, we see 
that our outer field still contains a fraction of sources belonging to the cluster or, more in 
general, associated to the recent star formation events in the Orion region (i.e. possible 
foreground members of the Orion OBI association). This is in agreement with what we 
found in the previous section discussing the color-magnitude diagrams, not by chance since 
the LF represents the projection on the vertical axis of the color-magnitude plot. 

The flattening of the substellar LFs in the inner region may be due to field star 
contamination. To address this possibility, we can use our outer field as a control field, 
knowing that it will probably provide an overestimate of the contribution of non-cluster 
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sources and therefore, after subtraction, a lower limit to the cluster LF. This for two main 
reasons: first, we have seen that our outer field still contains a fraction of cluster- related 
sources and, second, background sources should be more easily observed in the outer field 
rather than in the central region due to the larger background extinction caused by the 
underlying OMC-2, which has the ridge of highest density crossing the cluster core nearly 
in a north-sout h direction (see e . g. the extinction map obtai ned from the C 18 Q column 



densit y data of 



Goldsmith et al. 



Hillenbrand fc Carpenter 



(119971 ) presented in Fig 2 in 
(2000) ). Figure 23| (bottom, gray area) shows the LFs corrected in this way, having properly 
normalized the source densities to the same surface area. There is a clear reduction of 
source density in the substellar regime and the main peak becomes more prominent, but 
the substellar LFs seem to remain flat or show, at best, a very modest negative slope. 

Alternatively, it is possible to account for the higher background extinction beyond the 
central field, and therefore remove at least some of the systematic bias introduced by using 
the outer field, by shifting the outer LFs, for each fi l ter, b y magnitude values compatible 



with the interstellar reddening law of ICardelli et al. 



(119891 ) . This because with a proper, 



wavelength dependent, change of zero points one can use the outer field counts to reproduce 
the luminosity functions 
If we refer again to the 



the luminosity function s that would have been ob s erved with higher foreground extinction. 

~k3) 



Hillenbrand fc Carpenterl ((20001) extinction map, we see that the 
very central part of the OMC-2 is characterized by a visual extinction as high as Ay = 75 m . 
For a detaile d analysis one should account fo r the spatial variation of the extinction, as 
done e.g. by iHillenbrand fc Carpenterl (120001 ) . but for the moment let's simply assume 



the LFs of the outer region and put all sources under Aj ~ 10 mag, Ah — 6 mag, and 
Ak s — 4 mag, corresponding to Ay ~ 33 mag. It is immediate to see that shifting the 
monotonically increasing field LFs to the right by a few magnitudes and subtracting them 
from the central ones produces no visible change on the central ones. One can therefore 
safely conclude that the true LF of the inner region lies between the one derived above 
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subtracting the outer field source density (corresponding to a fully transparent background) 
and the one originally observed (corresponding to a fully opaque background), which means 
that the inner LFs overall remain quite flat in the brown-dwarfs regime (0. 075-0. 012M Q ). 
In a future paper we will address the implications of this result on the cluster Inititial 
Mass Function, evaluating the foreground extinction and the stellar parameters of the point 
sources on the basis of the rich set of multicolor photometry obtained for the HST Treasury 
Program on the Orion Nebula. 



7. Summary and conclusions 

We have presented a photometric survey of the Orion Nebula Cluster in the J-, H- and 
Ks-passbands carried out at the 4m telescope of Cerro Tololo. The survey, covering a field 
of about 30' x 40' centered about V southwest of the Trapezium, has been performed in 
parallel to visible observations of the same region made in La Silla (Da Rio et al. 2009). 
The two datasets constitute the first panchromatic survey covering simultaneously the 
Orion Nebula Cluster from the [/-band to the Kg-band. 

The final catalog, photometrically and astrometrically calibrated to the 2MASS system, 
contains 7759 sources, (includ ing 174 and 22 sourc es whose photometry has been taken 



directly from the 2MASS and 



Muench et al. 



(120021 ) catalogs, respectively). This represents 



the largest near-IR catalog of the Orion Nebula Cluster to date. Our sensitivity limits allow 
to detect objects of a few Jupiter masses under about Ay ~ 10 magnitudes of extinction. 

We present the color-color diagrams, color-magnitude diagrams and the luminosity 
functions for three regions centered on the Trapezium containing the same number of 
sources (excluding the M43 sub-cluster). Sources in the inner region typically show IR 
colors compatible with reddened T Tauri stars, whereas the outer fields are dominated by 
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field stars seen through an amount of extinction which decreases with the projected distance 
from the center. The color-magnitude diagrams allow to clearly distinguish between the 
main ONC population, spread across the full field, and background sources. After correction 
for completeness, the luminosity functions in the inner region remains nearly flat, with 
marginal contamination from background sources. 

We are grateful to the ESO director for allocating Director Discretionary Time with the 
ESO-MPIA 2.2m telescope in parallel with the ISPI time at Cerro Tololo. We are grateful 
to the ESO staff at La Silla for carrying out the WFI observations in service observing 
mode and to N. van der Bliek and the CTIO staff for support with the observations and 
data calibration. Observing time with ISPI was obtained within the framework of the HST 
Treasury Program GO10246. 

The reduced images presented in this paper can be downloaded from the High-Level 
Science Products archive hosted at the Multimission Archive at STScI at the address 

http: / / archive.stsci.edu/prepds/orion/ 



A. The database 



The following is a legend of the fields in the database. 



1. id: Running catalog number. The 174 sources taken f rom the 2MAS S cata log have a 



Muenchetal 



(120021 ) catalog 



progressive id starting from 10,000. Sources from the 
have their original entry number added to either 20,000 or 21,000, depending if the 
photometry was derived from FLWO or VLT observations (for sources having both 



measurements, the FLWO photometry has been adopted). 



2. ra: Source Right Ascension (2000.0) in decimal hours. Based on the Ks band value if 
the source is detected and not saturated, otherwise the H-band or the J-band, in the 
order. 

3. dec: Source Declination (2000.0) obtained following the same priorities of the RA 
above. 

4. raj: Right Ascension in the J band images. 

5. decj: Declination in the J band images. 

6. rah: Right Ascension in the Hband images. 

7. dech: Declination in the H band images. 

8. rak: Right Ascension in the K s band images. 

9. deck: Declination in the Ks band images. 

10. j: J band PSF magnitude, calibrated to the photometric 2MASS system. 

11. dj: Uncertainty in the above. 

12. h: Hband PSF magnitude, calibrated to the photometric 2MASS system. 

13. dh: Uncertainty in the above. 

14. k: K s band PSF magnitude, calibrated to the photometric 2MASS system. 

15. dk: Uncertainty in the above. 

16. j_ap: J band magnitude for a 5 pixel radius aperture, calibrated to the photometric 
2MASS system. 

17. dj_ap: Uncertainty in the above. 
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18. h_ap: H band magnitude for a 5 pixel radius aperture, calibrated to the photometric 
2MASS system. 

19. dh_ap: Uncertainty in the above. 

20. k_ap: Ks band magnitude for a 5 pixel radius aperture, calibrated to the photometric 
2MASS system. 

21. dk_ap: Uncertainty in the above. 

22. image: Image from which the measures are taken. Image numbers run from 1 to 11. 

23. exptime: Set to either 3 or 30, it refers to the exposure time in seconds of the 
individual frames building the image from which the photometry has been taken. 
Sources with exptime=3 are generally saturated in the 30 s exposures. 

24. status: Flag sources according to various criteria: 

— > 1 : source is point-like, detected during the first photometric extraction; 

— >• 2 : source is point-like, detected after PSF subtraction; 

— > 3 : source appears extened (galaxies, proplyds, diffuse objects in general); 

25. clone: Flag for multiple detections. Sources detected multiple times have the clone 
value equal to the id value of their first appearance in the database, multiplied by -1. 

Facilities: CTIO (ISPI) 
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g. 1. — Position of the eleven ISPI fields, superimposed to the positions of the sources 
tected in our survey. 
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Fig. 2. — Comparison between 2MASS magnitudes and ISPI calibrated magnitudes. 
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Fig. 3. — PSF photometric error as a function of magnitude for our point-like sources. The 
solid line indicates the locus of the average relation for the 2MASS sources falling in the 
ISPI field. 
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orcmin 

Fig. 4. — Projected regions identified in our surveyed area. The solid circles define the three 
radial areas containing an equal number of point-like sources, with radii of 6.7', 14.3' and 
27.2', respectively. The dashed circle indicates the exclusion field centered around NU Ori. 
The positions of 9 1 Ori-C, 9 2 Ori and NU Ori are also indicated, together with our original 
field positions. The axes show the angular offsets from 6 1 Ori-C (RA (2000.0)=05 /l 35 m 16.46 s , 
DEC (2000.0)=-05 o 23'23"18). 
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Fig. 5. — Completeness levels for the J (left), H (center) and K$ (right) bands derived from 
the artificial star experiment. Each plot shows the results for the 3 radial regions. 
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Fig. 6. — dashed lines: H and luminosity functions for the 



20001 ) (upper panels) and 



Muench et al. 



Hillenbrand fc Carpenter 



(]2002l ) (bottom panels) catalogs; gray area: lumi- 



nosity functions for the same filter and fields from the ISPI observations (point sources only); 
solid line: ISPI luminosity functions corrected for completeness. The error bars have been 
obtained by adding quadratically the statistical (poissonian) uncertainty on the measured 
counts to the completeness correction error derived from the Monte Carlo simulations. 
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Fig. 7. — Near-IR color composite mosaic of the Orion Nebula from our ISPI images. The 
image has been digitally enhanced to improve the visibility of the stellar sources in the 
central region. The RGB colors code the Ks, H and J bands, respectively. 
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Fig. 8. — Spatial distribution of the ONC stars for the inner (Trapezium) region. The 
caption at top of the figure illustrates how the size of the circles is proportional to the source 
brightness, whereas the colors are from blue to red according to their observed J-K$ index. 
White arrows are used to indicate the Trapezium stars, for orientation 
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Fig. 9. — Same as Figure [8], for Field 1. 
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Fig. 10. — Same as Figure for Field 2. 
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Fig. 11.- 



Same as Figured for Field 3. 
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Fig. 12.- 



Same as Figure for Field 4. 
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Fig. 13.- 



Same as Figured for Field 5. 
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Fig. 14- 



Same as Figured for Field 6. 
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Fig. 15- 



Same as Figure for Field 7. 
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Same as Figured for Field 8. 
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Fig. 17- 



Same as Figured for Field 9. 
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Fig. 18- 



Same as Figure for Field 10. 
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Fig. 19.- 



Same as Figured for Field 11. 
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Fig. 20. — a) top-left panel: Near-infrared color-color diagrams for selected point-like sources 
in our catalog detected in all 3 bands (see text for the selection criteria). Also shown are 
the main-seque nce (magenta solid line), the giant-branch (red solid line) sequence and the 



CTTs locus by 



Meyer et al. 



(119971 ) (green solid line). The dotted lines are parallel to the 
standard reddening vector, b) (top-right panel) - Hess diagrams of the color-color diagram 
relative to the inner region, c) bottom-left panel: Similar to panel b), for the median region. 
d) bottom-right panel: Similar to panel b) for the outer region. 
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Fig. 21. — Infrared color-magnitude diagrams for selected point sources in our catalog de- 
tected in all 3 bands, with lMyr is ochrones, assuming a distance d = 414 pc (DM = 8.01) . 



To derive these isochrones, we used (jChabrier et al. 



2000l ) combined to the (jSiess et al 



200(1 



model, as discussed in Section 16.21 The Ay =10 mag reddening vector is also shown. The 
diamond symbols represent the positions of a M = 0.012 M (deuterium burning limit) 
object, a M = 0.075 M (hydrogen burning limit) object and a M = 1 M star from bottom 
to top respectively. 
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Fig. 22. — Hess format for the J,J-H (left panels) and Ks,H-K$ (right panels) color- 
magnitude diagrams. Line styles and symbols are the same as in Figure [2TJ From top 
to bottom, the diagrams for the inner (r < 6.7'), medium (6.7' < r < 14.3') and outer 
(14.3' < r < 27.2') regions are shown. 
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Fig. 23. — Top panels: J-band (left), ff-band (center) and K^-band (right) luminosity 
functions for the ONC, corrected for completeness. The solid line is relative to the inner 
region, the dashed line to the intermediate region and the dotted-dashed lines to the outer 
region. Bottom panels: for each bandpass, luminosity function of the inner region (solid 
line), of the outer region shifted by Ay = 33 mag of extinction (dot-dashed line), and 
difference of the two (gray area). In each plot, we indicate with vertical lines the magnitudes 
corresponding to stars 1 Myr old at 414 pc with 1 M , 0.075 M Q and 0.012 M respectively, 
without interstellar reddening. The vertical dashed line stands for the 80% sensitivity level 
in the outer region. 
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Table 1. ISPI Broad Band Filters 



Filter 


Central wavelength 


AA(80%) 




Om) 


H 


J 


1.250 


1.176 - 1.322 


H 


1.635 


1.500 - 1.770 


K S 


2.150 


1.992 - 2.300 
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Table 2. Log of the ISPI observations 



Field Night RA (2000.0) DEC (2000.0) Airmass Comment 



1 


A 


05 


:35:38.37 


-5:12:50.7 


1.144 - 


1 


.291 




2 


A 


05 


:34:59.63 


-5:12:57.6 


1.123 - 


1 


.139 




3 


A 


05 


:34:18.52 


-5:12:59.6 


1.123 - 


1 


.139 


no J-band 30s exposure 


4 


A 


05 


:35:37.05 


-5:22:38.6 


1.309 - 


1 


.158 






B 


05 


:35:36.96 


-5:22:50.6 


1.309 - 


1 


.158 


field 12 in the database 


5 


A 


05 


.34:59.08 


-5:22:28.3 


1.319 - 


1 


.938 






B 


05 


:34:59.21 


-5:22:55.0 


1.225 - 


1 


.498 


field 13 in the database 


6 


A 


05 


:34:17.98 


-5:22:28.3 


1.319 - 


1 


.938 






B 


05 


.34:18.01 


-5:22:55.7 


1.225 - 


1 


.498 


field 14 in the database 


7 


B 


05 


:35:38.16 


-5:33:03.8 


1.218 - 


1 


.410 




8 


B 


05 


:34:59.82 


-5:33:18.6 


1.104 - 


1 


.211 




9 


B 


05 


:34:18.50 


-5:33:21.4 


1.104 - 


1 


.211 




10 


B 


05 


.35:37.48 


-5:43:06.9 


1.110 - 


1 


.215 




11 


B 


05 


:34:56.28 


-5:43:09.2 


1.110 - 


1 


.215 
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Table 3. Best estimate for calibration coefficients. 



Field 


Z H 






ff e.JK 


1 


1.09 


0.01 


-0.010 


0.008 


2 


1.07 


0.01 


0.001 


0.007 


3 


1.08 


0.04 


0.00 


0.02 


4 


1.13 


0.01 


-0.022 


0.006 


5 


1.11 


0.01 


-0.023 


0.007 


6 


1.03 


0.02 


-0.05 


0.02 


7 


0.93 


0.02 


0.00 


0.01 


8 


0.93 


0.01 


-0.023 


0.008 


9 


0.95 


0.01 


-0.02 


0.01 


10 


1.09 


0.02 


-0.06 


0.01 


11 


1.07 


0.02 


-0.02 


0.01 



Note. — Parameters for field 3 are com- 
puted using non-normalized J-band magni- 
tudes. 
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Table 4. Source Summary 



Source Number of stars 



ISPI detected in J-, H- and K^-bands 4851 

ISPI detected in J- and if-bands, not in K^-band 92 

ISPI detected in H- and K^-bands, not in J-band 1096 

ISPI detected in J- and .FCs-bands, not in ff-band 85 

ISPI detected only in J-band 175 

ISPI detected only in ff-band 80 

ISPI detected only in Kg-band 251 

ISPI diffuse 933 

ISPI Total 7563 

Extracted from the 2MASS catalog 174 

Extracted from the Muench et al. (2002) catalog 22 

Grand Total 7759 

ISPI detected in 30 s images 6147 

ISPI detected in 3 s images 1416 

ISPI first run detection 6305 

ISPI hidden companions 325 

ISPI diffuse 933 
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Table 5. ISPI Photometry and Astrometry of the Orion Nebula Cluster 



ID 


RA 


DEC 


RAJ 


DECJ 


RAH 


DECH 


RAK 


DECK 


.1 


D.I 


H 



1 


5. 


.5891618 


-5 


.2200644 


0. 


.0000000 


0. 


0000000 





.0000000 


0. 


.0000000 


5. 


5891618 


-5. 


2200644 


0.00000 


0.00000 


0.001 


2 


5. 


.5953268 


-5 


.2236077 


5. 


.5953268 


-5. 


.2236077 


5 


.5953268 


-5. 


.2236077 


5. 


5953268 


-5. 


2236077 


0.00000 


0.00000 


0.001 


3 


5. 


.5892962 


-5 


.2290140 


5. 


.5892962 


-5. 


.2290140 


5 


.5892962 


-5. 


.2290140 


5. 


5892962 


-5. 


2290140 


0.00000 


0.00000 


0.001 


1 


5. 


.5957626 


-5 


.2373337 


5. 


.5957626 


-5. 


.2373337 


5 


.5957626 


-5. 


.2373337 


5. 


.5957626 


-5. 


.2373337 


0.00000 


0.00000 


0.001 


5 


5. 


.5977591 


-5 


.2379900 


5. 


.5977591 


-5. 


.2379900 


5 


.5977591 


-5. 


.2379900 


5. 


5977591 


-5. 


2379900 


0.00000 


0.00000 


0.001 


6 


5. 


.5977757 


-5 


.2385170 


5. 


.5977757 


-5. 


.2385170 


5 


.5977757 


-5. 


.2385170 


5. 


5977757 


-5. 


2385170 


0.00000 


0.00000 


0.001 



50 


5. 


.5892075 


-5 


.3059220 


5. 


.5892064 


-5. 


.3059297 


5 


.5892100 


-5. 


.3059316 


5. 


5892075 


-5 


.3059220 


13.0184 


0.0607901 


11.9: 


51 


5. 


.5884842 


-5 


.3057194 


5. 


.5884816 


-5. 


.3056698 


5 


.5884806 


-5. 


.3056545 


5. 


5884842 


-5 


.3057194 


12.7213 


0.0437424 


11.9! 


52 


5. 


.5919210 


-5 


.3049874 


5. 


.5919174 


-5. 


.3049865 


5 


.5919189 


-5. 


.3049850 


5. 


.5919210 


-5 


.3049874 


11.4727 


0.0458089 


10.71 


53 


5. 


.5919566 


-5 


.3021426 


5. 


.5919551 


-5. 


.3021808 


5 


.5919540 


-5. 


.3021369 


5. 


5919566 


-5 


.3021426 


17.4306 


0.218837 


16.6: 



3597 5.5816976 -5.4131470 5.5817022 -5.4131479 5.5816945 -5.4131670 5.5816976 -5.4131470 16.7388 0.112858 16.01 

3678 5.5823130 -5.4152808 0.0000000 0.0000000 5.5823110 -5.4152832 5.5823130 -5.4152808 0.00000 0.00000 18.0: 

4109 5.5791845 -5.3399510 5.5791916 -5.3398418 5.5791794 -5.3398919 5.5791845 -5.3399510 17.3715 0.134945 16.3! 

82996 5.5871119 -5.5159988 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 10.1820 0.0230000 9.26: 

82620 5.5860662 -5.4647789 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 7.74300 0.0240000 7.631 

81916 5.5811772 -5.5523648 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 8.45700 0.0320000 7.951 



81767 5.5797176 -5.5707202 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 0.0000000 7.22100 0.0190000 6.96- 



Note. - 



Tabic [5] is published in its entirety in the electronic edition. A portion is shown here for guidance regarding its form and content. See 



